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Basic properties of the nuclear tensor mean fields are reviewed, and their role in 
changing the shell structure and masses of nuclei is analyzed within the spherical 
Hartree-Fock-Bogolyubov approach. 



1. Introduction 



Modern formulation of the nuclear mean-field theory is based on the en- 
q ■ ergy density formalism, 1 which has been over the years developed for elec- 

^ . tronic systems. According to the formal Hohenberg-Kohn 2 and Kohn-Sham 

theorems, 3 exact ground-state energies of many-fermion systems can be ob- 
tained by minimizing certain exact functional of one-body density. These 
theorems do not provide any method to construct the exact functional in 
Vh ' a systematic way; nevertheless one can build phenomenological functionals 

and test their performance against experimental data. Such an approach is 
also consistent with the ideas of the effective field theory, whereupon prop- 
erties of composite objects at low energies can be described by Lagrangians 
which include high-energy dynamics in the form of the appropriate series 
of contact terms. 

Within the energy density formalism, one treats the nucleus as a single 
composite object described by a set of one-body densities. At low energies, 
when the densities are varying slowly in the nuclear interior and then go 
smoothly to zero at the nuclear surface, one can consider only local densities 
that are built of the one-body density matrix and its derivatives up to the 
second-order. Systematics construction of the most general energy density 
functional (EDF) consistent with symmetries is then possible, 4 and gives 
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a generalization of the extremely successful approach based on the Skyrmc 
effective interaction. 5 

The origins of the spin-orbit (SO) splitting in nuclei can be attributed to 
the bare two-body SO and tensor interactions, 6 which contribute differently 
to the spin-saturated (SS) and spin-unsaturated (SUS) nuclei. An alterna- 
tive explanation is also sought in the relativistic mean-field theories with 
meson couplings 7 , where no distinction between the SS and SUS systems 
is obtained. 

The Skyrme interaction was introduced into nuclear physics more then 
30 years ago, 8 and shortly after it was supplemented by tensor forces. 9,10 
However, after these ground-breaking studies, in most of the subsequent 
applications the tensor forces were not taken into account. Moreover, within 
many Skyrme-force parameterizations constructed to date, the tensor terms 
in the EDF that were coming from the central force, were quite arbitrarily 
set equal to zero. 5 

In the present paper, I discuss the form of the tensor terms in the EDF 
and their influence on the single-particle energies and nuclear masses. In 
Sec. 2, I recall the recent experimental evidence on the changes of shell 
structure in neutron-rich Z«20 nuclei. 11,12 This is only one of several such 
examples recently identified in light nuclei and interpreted within the shell- 
model by introducing tensor interactions. 13,14 Properties of the tensor terms 
in the EDF are discussed in Sec. 3, and in Sec. 4, I present results of 
calculations for single-particle energies and masses obtained with tensor 
terms included in the mean-field approach. 



2. Shell structure of neutron-rich Z~20 nuclei 

In a series of recent experiments performed at the Argonne National Lab- 
oratory with Gammasphere 11 and National Superconducting Cyclotron 
Laboratory, 12 properties of low- lying collective states of even-even neutron- 
rich titanium isotopes have been measured. As illustrated in Fig. 1, the data 
reveal the presence of a closed A=32 subshcll, in addition to the standard 
A^=28 shell present in heavier elements. Indeed, both 50 Ti and 54 Ti show 
the increased 2 + energies, and decreased BE2 values, as compared to their 
neighbours. 

In order to explain such a change of the shell structure, the single- 
particle neutron viz/2 orbital must be shifted up, which leaves a gap be- 
tween the spin-orbit-split ^P3/2 and ^Pi/2 orbitals, and creates a subshcll 
closure for four particles occupying vp^/2- The shell- model calculations, 15 
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Figure 1. Energies of the first excited 2+ states (squares, left scale) and reduced tran- 
sition probabilities BE2(0+^2+) (circles, right scale), measured in neutron-rich Ti 
isotopes. 11,12 

performed for the single-particle orbitals shifted in this way, confirm the 
pattern shown in Fig. 1. The origins of the shift are attributed to the de- 
creased monopole interaction energy between the proton irf>i/2 and neutron 
^5/2 orbitals, which occurs when protons are removed from 7rf7/ 2 - The 
source of such a monopole interaction is in turn attributed to the shell- 
model tensor interaction between these orbitals. 

Positions of single-particle levels can be best studied within the mean- 
field approximation, in which they are basic dynamic characteristics of the 
system, resulting from the two-body interactions being averaged with par- 
ticle densities of occupied states. Therefore, in this paper I evaluate the 
single-particle energies by applying the mean-field methods to tensor inter- 
actions. 



3. Tensor densities in the energy density functional 

3.1. Spin-orbit and tensor forces 

Momentum-dependent two-body SO 8 and tensor 9 ' 10 interactions have the 
form 



iW S ■ 


k x k 
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k ■ S • k +k ■ S • k 
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where the vector and tensor spin operators read 
S = <J\ 



S' J 



<T2, 



(5 U <Tl • <T 2 . 



(2) 



When averaged with one-body density matrices, these interactions con- 
tribute to the following terms in the EDF (see Refs. 16 ' 4 for derivations), 

~Hso = jWo 3J - 



r dpi 

Jl dr 



(3) 



where the conservation of time-reversal and spherical symmetries was as- 
sumed. Here, pt and Jt are the neutron, proton, isoscalar, and isovector 
particle and SO densities 8, 16,4 for t=n, p, 0, and 1, respectively. 

Apart from the contribution of the SO energy density to the central 
potential, variation of the SO and tensor terms with respect to the densities 
yields the one-body SO potential for neutrons (t=n) and protons (t=p), 



W t so 



It 



W 



dr 



+ l£)+l((te+t )Jo-(te~t )J t 



L-S. 



(4) 



Hence, it is clear that the only effect of including the tensor interaction 
is a modification of the SO splitting of the single-particle levels, and that, 
from the point of view of one-body properties, tensor interactions act very 
similarly to the two-body SO interactions. However, the latter ones induce 
the SO splitting that is only weakly depending on the shell filling. This is 
so because the corresponding form-factor in Eq. (4) is given by the radial 
derivatives of densities. On the other hand, the SO splitting induced by the 
tensor forces depends strongly on the shell filling, because its form-factor is 
given by the SO densities J(r). Indeed, when only one of the SO partners 
is occupied (SUS system), the SO density is large, and when both partners 
are occupied (SS system), the SO density is small, see Sec. 4 for numerical 
examples. 



3.2. Spin-orbit and tensor energy densities 

Within the energy-density approach, one does not relate the EDF to an 
average of the two-body force, but one postulates the EDF based on sym- 
metry conditions only. Then, the most general EDF, depending on the 
spin-current densities, reads 4 

n S o = Et=o,iC7 J ptV-J t , 



n T = E*=o,i ( c t°J" + Ct J1 J 2 t + Ct J2 ±t 



(5) 
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where the spherical-symmetry condition has been released. The standard 
pseudoscalar J_ t , vector J t , and pscudotensor J t parts of the spin-current 
density, 16 

(6) 
r=r' 





T 1 

•Jabt — 2 i 


(V-VVw(r,r ; ) 


defined as 






d-t = l^a=x,y,z ■'ooti 




J at = 


Z--ib,c=x,y,z e abcJbcti 



(7) 

—abt — 2 a °t ' 2 b a t ~^lLt^ab- 

The tensor energy density TLt now depends on six coupling constants, 
C/°, C/ 1 , and C/ 2 , for £=0,1, and not on two coupling constants, t e and 
t , as in Eq. 3. Similarly, the SO energy density Tiso depends now on 
two coupling constants C t VJ for £=0,1, and not on one, Wo, (the latter 
generalization has been introduced and studied in Rcf. ). 

From the symmetry conditions imposed by the spherical, axial, and 
reflection symmetries one obtains 18 that: 

1° The pseudoscalar densities J_ t vanish unless the axial or reflection 

symmetries are broken. 
2° The pseudotensor densities J t vanish unless the spherical symmetry 

is broken. 

Finally, the gauge-invariance symmetry conditions 4 require that there are 
only two gauge-invariant combinations of the pseudoscalar, vector, and 
pseudotensor terms, namely, 

"t = 3— t + 2 J t + It ' . . 

rF _ 2 t2 _ 1 j2 , 1 ,2 ^ ' 

^t — 3^-t 4 J t T" 2 ±t ■ 

In such a case, only four out of the six tensor coupling constants are linearly 
independent, i.e., 

C t = g A t + g-Bt, 

C/ 1 = \A t - \B t , (9) 

C/ 2 = A t + \B t . 

On the other hand, the averaging of the tensor forces (1) implies that only 
two out of the six tensor coupling constants remain linearly independent, 
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i.e., 



Bq — — 3Ao — — g (t e + 3t ) 
B\ = —3Ai = -° (t e — to) ■ 



(10) 
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Figure 2. Radial components J n (r) of the neutron vector SO densities J n , calculated 
for the Z=28 isotopes with N =28-38 (left panel) and JV=40-50 (right panel). 



In order to illustrate the influence of tensor densities on single-particle 
and global nuclear properties, in Fig. 2 are shown the radial components, 
Jn{r)=Jn • r / r i °f the neutron vector SO densities J n , calculated for the 
nickel isotopes between N—28 and 50. Calculations have been performed 
for the Skyrme SLy4 interaction, 19 by using the Hartree-Fock-Bogolyubov 
(HFB) method with the spherical symmetry assumed. 20 

One can see that the SO densities are mostly positive and peaked near 
the surface. At 7V=28, the SO density is large and its major part comes 
from the occupied ^7/2 orbital (SUS system). By adding neutrons in the 
shell above N—28, this part is gradually cancelled by an increasing in mag- 
nitude, negative contribution from the SO partner 2^5/2- At the same time 
contributions from the ^3/2 and z/pi/2 orbitals appear. When both pairs 
of the SO partners are occupied around 7V=38, and when the z^gg/2 orbital 
is still empty (SS system), the SO density is rather small. Beyond iV=40, 
it increases again until the z^gg/2 orbital becomes fully occupied at A=50. 
Note the shift of the SO densities to larger distances, which occurs at the 
point of the switch-over between the dominating 1/ and \g contributions. 
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A similar pattern of varying SO densities is valid for all shells. For SS 
systems, one obtains small SO densities, while for SUS systems, the SO 
densities are large. Therefore, the SO densities are small at magic shells 
N, Z=2, 8, and 20, and large at magic shells N, Z=28, 50, 82, and 126. 

Since the j > partners are always occupied first, the SO densities are 
mostly positive. Note that the derivatives of particle densities are mostly 
negative, and also peaked at the surface; therefore, for positive coupling 
constants, the SO and tensor forces split the SO partners in opposite direc- 
tions, cf. Eq. (4). 

4.2. Single-particle levels 
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Figure 3. Single-particle properties of neutron levels in N=32 isotoncs, calculated 
within the HFB method with the SLy4 Skyrmc interaction. Left panels correspond to 
the standard SLy4 paramctrization with no tensor terms (i e =0) while the right panels 
correspond to the tensor-even interaction included (4 e =200McVfm 5 ). Single- particle 
energies in the sdpf shell (top panels) arc shown along with their spin-orbit-splitting 
energies (middle panels) and centroid energies (bottom panels). 



In Fig. 3 are shown properties of neutron single-particle levels calculated 
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for the chain of the 7V=32 isotones (these levels are relevant for the changes 
of the shell structure discussed in Sec. 2). Single-particle energies e 7 ^ 3 (top 
panels) were calculated as the canonical energies 21 of the HFB method. In 
order to better visualize the dependence of the single-particle energies on 
the proton number, in the middle and bottom panels are shown the SO 
splittings and centroids, respectively, of the SO partners, defined as 



£ ccnt 2 \ £ Af 



1 ( ntj < . nlj~ 



(11) 



Without the tensor terms (left panels of Fig. 3), the neutron single- 
particle energies vary smoothly with the proton numbers. Apart from a 
gradual increase with decreasing Z, one observes two clear type of changes 
in the shell structure. First, in each shell the centroids of levels with dif- 
ferent values of i become degenerate towards the neutron drip line. This 
effect is related to the increase of the surface diffuseness of particle distribu- 
tions, which renders the shell structure of very neutron-rich nuclei similar 
to that of a harmonic oscillator. 22 Second, near the neutron drip line the 
SO splitting of the weakly bound p orbitals becomes smaller, because such 
orbitals start to decouple from the SO potential due to their increasing 
spatial dimensions. 23 

In the right panels of Fig. 3 are shown the analogous results ob- 
tained with the tensor-even interaction Vr e , Eq. (1), taken into account 
for £ e =200MeVfm 5 . (This particular value of the coupling constant was 
not optimized in any sense, and it is used here only to illustrate some general 
trends.) The use of the tensor-even interaction (T=0, S=l neutron-proton 
channel) corresponds to the energy density that depends on the product of 
neutron and proton SO densities, cf. Eq. (3). Therefore, the effect of the 
tensor term vanishes at the closed proton shell Z=20 (SS system). For Z 
higher (lower) than 20, the effect of the tensor term increases as a result 
of increasing contributions to the proton SO density coming from the Tri^/2 
(7^3/2) orbitals. This is clearly visible in the middle right panel of Fig. 3, 
where the SO splitting decreases on both sides of 52 Ca. This is so because, 
for positive values of the coupling constant t e , the effect of the tensor force 
partly cancels that of the standard SO force. As a result, the 7rfs/2 orbital 
is in 60 Ni much closer to the p orbitals than it is in 54 Ti; the shift which is 
compatible with the changes of the shell structure discussed in Sec. 2. 

A detailed reproduction of the level positions is not the goal of the 
present study. The coupling constants of the tensor terms have to be ad- 
justed together with other parameters of the EDF, by considering not only 
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this particular region of nuclei, and not only this particular set of observ- 
ables. Indeed, the tensor terms included in the EDF will influence many 
different global nuclear properties throughout the mass chart, and a global 
analysis is therefore necessary. Before this is done, in the next section, the 
impact of the tensor interaction on nuclear binding energies is studied in a 
preliminary way. 

4.3. Binding energies 

When the tensor terms (3) are added to the EDF, the binding energies are 
affected through self-consistent changes of all the terms in the EDF. How- 
ever, qualitatively, the effects of tensor terms on the ground-state energies 
can be illustrated by integrals of products of the SO densities that appear 
in Eq. (3). In Fig. 4, values of such integrals are shown for the neutron 
SO densities squared, J„(r), calculated at magic proton numbers in func- 
tion of the neutron numbers. Comparison of results obtained without (left 
panel) and with (right panel) tensor-even interaction included, shows that 
the effect of the tensor term can, in the first approximation, be treated per- 
turbatively. Due to the fact that the proton SO densities depend weakly 
on the neutron numbers, for .2=28, 50, and 82 the integrals of products 
Jn( r )Jp{f) show similar a behaviour to those of J^(r), while they are small 
for Z=8 and 20. 

From the results shown in Fig. 4, it is clear that, for positive cou- 
pling constants, the tensor terms will give characteristic contributions to 
the ground-state energies of heavy nuclei. These contributions will have 
a form of inverted arches, spanned between the neutron magic numbers. 
This feature is conspicuously reminiscent of differences between the theo- 
retical and experimental ground-state energies obtained with tensor terms 
not included. 24 ' 25 ' 26 It is therefore quite plausible that by including the 
tensor terms one may be able to remove a major part of the discrepancy 
between the previously calculated nuclear masses and experiment. 

5. Conclusions 

During the past thirty odd years, when the nuclear self-consistent mean- 
field methods based on effective interactions were developed and imple- 
mented, the tensor interactions have been largely ignored. On the other 
hand, recent experimental studies and shell-model analyses indicate that 
these interactions may play an important role in several regions of nuclear 
chart. A unified picture of the role played by tensor interactions through- 
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Figure 4. Integrals of the neutron SO densities squared that define contributions of the 
tensor terms to total binding energies. Left panels correspond to the standard SLy4 
paramctrization with no tensor terms (£ e =0), while the right panels correspond to the 
tensor-even interaction included (£ e =200McVfm 5 ). 

out the mass table is still missing, and is very much needed. In the present 
paper, I have reviewed basic properties of the tensor mean fields, and I have 
illustrated their role in changing the shell structure and masses of nuclei. 
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